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Distorsiones  aleatorias  en  la  adquisición  de  la  imagen 


Introducción 


La  imagen  adquirida  por  un  sensor  a  bordo  de  un 
satélite  presenta  diferencias  con  respecto  a  la  imagen 
real  tanto  en  la  situación  de  sus  puntos  como  en  el 
nivel  de  brillo  de  cada  uno.  En  consecuencia,  la  geo¬ 
metría  de  la  imagen  no  refleja  con  fidelidad  absoluta 
la  de  la  superficie,  ni  tampoco  pueden  deducirse  sus 
propiedades  físicas  a  partir  de  la  lectura  directa  de 
los  niveles  digitales.  Es  necesario  aplicar  unas  deter¬ 
minadas  operaciones  de  corrección,  que  están  enca¬ 
minadas  a  minimizar  estas  alteraciones  en  la  medida 
de  lo  posible.  Algunas  de  estas  tareas  son  aplicadas 


Distorsiones 


Las  imágenes  de  satélite  no  constituyen  una  carto¬ 
grafía  debido  a  factores  como  la  propia  naturaleza  de 
la  proyección,  las  distorsiones  sistemáticas  y  las  alea¬ 
torias  que  sufren  las  imágenes  como  consecuencia  de 
la  geometría  de  la  adquisición,  del  movimiento  del 
satélite  y  de  la  óptica  y  la  mecánica  del  instrumento 
de  adquisición.  Las  distorsiones  que  sufre  una  ima¬ 
gen  pueden  ser  de  tipo  geométrico  o  radiométrico. 


Distorsiones  geométricas 

Las  fuentes  que  dan  origen  a  las  distorsiones  geomé¬ 
tricas  se  enumeran  a  continuación. 


por  las  estaciones  de  recepción  o  los  distribuidores 
de  las  imágenes,  que  ofrecen  distintos  niveles  de  co¬ 
rrección.  El  estudio  detallado  de  la  imagen,  los  estu¬ 
dios  multitemporales  o  la  fusión  con  otras  fuentes 
cartográficas  requieren  correcciones  precisas,  que 
normalmente  son  realizadas  por  el  analista. 

Este  capítulo  aborda  las  operaciones  más  habituales 
en  el  pretratamiento  de  las  imágenes  de  satélite  antes 
de  proceder  al  estudio  en  sí  de  la  información  que 
proporcionan. 


Oscilación  de  la  plataforma 

Si  bien  el  movimiento  orbital  del  satélite  es  práctica¬ 
mente  uniforme,  y  por  tanto,  sus  aceleraciones  linea¬ 
les  son  nulas  en  cualquier  dirección  salvo  las  deriva¬ 
das  de  la  excentricidad  de  su  órbita,  en  realidad 
puede  sufrir  pequeñas  variaciones  en  altura  orbital, 
velocidad  u  orientación  de  cualquiera  de  sus  tres 
ejes,  lo  que  da  lugar  a  los  efectos  de  aleteo  (pitch), 
que  es  el  giro  según  el  eje  de  la  trayectoria,  de  cabe¬ 
ceo  (ro//),  o  giro  según  un  eje  perpendicular  al  pri¬ 
mero  y  paralelo  a  la  superficie  terrestre,  y  ladeo 
iyaw),  giro  según  el  eje  perpendicular  a  ella. 

Rotación  terrestre 

La  traza  de  un  satélite  polar  no  está  situada  a  lo  largo 
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de  un  mismo  meridiano  terrestre,  debido  a  que  mien¬ 
tras  el  satélite  describe  su  órbita,  la  Tierra  rota  alre¬ 
dedor  de  su  eje.  Cuando  el  satélite,  tras  haber  cru¬ 
zado  el  polo  terrestre,  ha  alcanzado  el  ecuador,  el 
planeta  ha  girado  y  el  punto  ecuatorial  que  pertenece 
al  meridiano  contenido  en  el  plano  orbital  cuando  el 
satélite  estaba  sobre  la  latitud  90-  se  ha  desplazado 
hacia  el  este. 

Los  satélites  adquieren  lineas  completas  de  imagen 
sucesivamente  y  sin  solución  de  continuidad  mien¬ 
tras  se  desplazan  en  su  órbita  alrededor  de  la  Tierra. 
La  escena  formada  por  un  cierto  paquete  de  lineas  se 
forma  en  el  segmento  terrestre  de  la  misión  para  co¬ 
mercializar  la  información.  Una  escena  contiene  un 
determinado  número  de  líneas  de  imagen,  que  son 
perpendiculares  a  la  traza  del  satélite.  Si  la  imagen 
corresponde  a  una  traza  descendente,  las  líneas  infe¬ 
riores  corresponderán  a  zonas  geográficas  situadas 
más  al  oeste  que  las  superiores,  que  fueron  tomadas 
unos  instantes  antes.  Lo  contrario  ocurre  si  la  traza 
es  ascendente.  Esto  produce  un  efecto  típico  de  dis¬ 
torsión  en  todas  las  imágenes  de  sensores  polares  o 
casi-polares  (no  afecta  a  los  geoestacionarios):  si  la 
imagen  no  se  corrige,  tendrá  forma  rectangular  y  las 
zonas  inferiores  estarán  desplazadas  hacia  el  este 
(traza  descendente).  Tras  la  corrección  de  la  distor¬ 
sión  debida  a  la  rotación  terrestre,  la  imagen  adqui¬ 
rirá  una  forma  romboidal  con  las  líneas  inferiores 
desplazadas  hacia  el  oeste,  representando  zonas, 
ahora  sí,  en  su  posición  geográfica  correcta. 

Tiempo  de  barrido 

Si  el  sensor  es  del  tipo  escáner,  cada  línea  o  grupo 
de  ellas  es  adquirida  por  un  solo  detector  o  un  pe¬ 
queño  número  de  éstos  mientras  un  espejo  explora 
toda  la  línea.  En  el  barrido  se  invierte  un  tiempo  de¬ 
terminado,  durante  el  cual,  el  satélite  evoluciona  en 


su  propia  órbita  y  deja  de  encontrarse  sobre  la  verti¬ 
cal  del  mismo  punto  en  el  que  comenzó  el  barrido.  El 
final  de  la  línea  corresponde  a  puntos  del  terreno  más 
hacia  el  sur  que  los  del  inicio  (en  traza  descendente). 
Es  otra  distorsión  que  debe  corregirse  para  hacer  útil 
la  imagen  desde  el  punto  de  vista  cartográfico. 

Distorsión  panorámica 

Dado  que  el  campo  de  visión  instantáneo  es  una  ca¬ 
racterística  propia  del  sensor  y  además  es  invariable, 
el  tamaño  de  las  celdas  en  puntos  alejados  de  la  traza 
del  satélite  es  mayor  que  el  correspondiente  al  nadir 
del  instrumento,  que  coincide  con  su  resolución  es¬ 
pacial  nominal.  Realmente,  ni  siquiera  el  IFOV  es  cua¬ 
drado  ni  rectangular,  sino  circular  o  elíptico.  La  for¬ 
mar  la  imagen  se  traduce  el  nivel  de  radiancia  regis¬ 
trado  en  cada  área  explorada  instantáneamente  a  un 
nivel  digital  conferido  a  una  celda  cuadrada. 

Curvatura  de  la  Tierra. 

Además  del  efecto  anterior,  hay  que  considerar  que 
la  Tierra  no  es  plana  y  que,  por  lo  tanto,  el  tamaño 
de  las  celdas  alejadas  de  la  traza  crece  considerable¬ 
mente.  Este  efecto  es  particularmente  apreciable  en 
satélites  meteorológicos,  de  gran  campo  de  visión,  so¬ 
bre  todo  en  los  geoestacionarios. 


Distorsiones  ra(diométricas. 

En  el  otro  grupo  de  distorsiones  se  encuentran  las  de 
tipo  radiométrico,  que  no  afectan  a  las  posiciones  de 
los  puntos,  sino  a  su  nivel  digital.  Normalmente  se 
trata  de  pérdida  de  celdas  aisladas  o  de  líneas  com¬ 
pletas,  que  aparecen  en  la  imagen  en  negro  o  en 
blanco,  o  bien  en  niveles  de  gris  anormalmente  altos 
o  bajos  con  respecto  a  su  entorno.  El  efecto  de  líneas 
alteradas  puede  repetirse  regularmente  a  lo  largo  de 


movimiento 
de  io  Tierra 


despiozo  miento 
relativo  de  io 
piotoforma 


sentido  del  barrido 


P 


Distorsión  debida  a  la  rotación  terrestre  (arriba),  al 
tiempo  de  barrido  (en  el  centro)  y  panorámica  (abajo). 


Carlos  Pinilla 


Operadores 


08-3 


Distorsión  debida  da  la  curvatura  terrestre. 


la  imagen  en  sensores  de  barrido  cuando  se  descali¬ 
bra  algún  detector,  causando  el  efecto  conocido 
como  bandeado.  En  los  sensores  de  empuje,  la  des¬ 
calibración  de  los  detectores  no  se  presenta  en  las 
líneas  sino  en  las  columnas  de  la  imagen  y  además 


no  suelen  seguir  ningún  patrón  repetitivo.  Final¬ 
mente,  cabe  también  hablar  de  una  distorsión  radio- 
métrica  causada  por  el  efecto  atmosférico,  ya  que 
tanto  la  dispersión  como  la  absorción,  aerosólica  o 
molecular,  alteran  los  niveles  digitales,  produciendo 
en  la  imagen  niveles  de  brillo  que  no  son  representa¬ 
tivos  de  la  radiancia  de  la  superficie. 


Operadores 


Un  operador  es  una  función  matemática  que  se  aplica 
sobre  la  imagen,  entendida  como  una  matriz  numé¬ 
rica.  Un  operador  T  es  una  aplicación  entre  el  espacio 
inicial,  por  ejemplo,  una  imagen  original  z=f(x,y)  y  el 
espacio  final,  constituido  por  la  imagen  transformada 
z'=f(x',y'): 

z  =  f{x,y)  ^  >z=f'{x',y') 

Esta  transformación  se  puede  desglosar  en  sus  tres 
componentes,  donde  la  variables  x  e  y  y  sus  transfor¬ 
madas  x'  e  y'  hacen  referencia  a  las  dimensiones  pla¬ 
nimétricas  de  la  imagen  y  la  variable  z  y  su  transfor¬ 
mada  z '  se  refieren  al  nivel  digital  de  la  celda: 

x'  =  t,(x,y,z) 
y' =  ty{x,  y,z) 
z=t,{x,y,z) 

Las  dos  primeras  constituyen  el  grupo  de  transforma¬ 
ciones  geométricas,  que  alteran  la  geometría  de  la 
imagen  sin  modificar  los  niveles  digitales  —o  al  menos 
no  de  forma  significativa,  como  se  verá  más  ade¬ 


lante-  y  la  tercera  las  transformaciones  radiométri- 
cas,  que  modifican  el  brillo  sin  intervenir  en  la  geo¬ 
metría.  Habitualmente,  ambos  tipos  de  transformacio¬ 
nes  se  abordan  por  separado.  Las  transformaciones 
también  pueden  escribirse: 

/i(x',y')  =  ri[/(x,y)] 
f2{z)  =  T^[f{z)\ 

Siendo  T\  y  Ti  lineales  e  invariables  al  desplaza¬ 
miento,  esto  es,  operan  del  mismo  modo  sobre 
cualquier  celda  de  la  imagen,  manteniéndose  la  fun¬ 
ción  igual  a  sí  misma: 

r[/(x-Xo,y-yo)]-/(x-Xo,y-yo)  = 

=  r[/(x,y)]-/(x,y) 

Según  el  número  de  celdas  involucradas  en  la  trans¬ 
formación,  los  operadores  pueden  ser: 
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Puntuales:  cuando  cada  elemento  de  la  ima¬ 
gen  final  es  función  de  un  único  elemento 
de  la  imagen  inicial. 

Locales:  cuando  el  nivel  digital  de  cada  ele¬ 
mento  de  la  imagen  final  es  función  de  los 
valores  de  las  celdas  de  un  determinado  en¬ 
torno  de  cada  elemento  de  la  imagen  inicial. 
Globales:  cuando  la  imagen  final  en  su  totali¬ 
dad  es  función  de  toda  la  imagen  inicial. 


Entre  los  primeros  se  encuentran  unas  de  las  opera¬ 
ciones  más  habituales  en  el  tratamiento  digital  de 
imágenes,  las  transformaciones  radiométricas  y  geo¬ 
métricas.  Las  transformaciones  radiométricas  modifi¬ 
can  la  amplitud  de  la  señal  (brillo  o  nivel  digital), 
manteniendo  invariables  las  coordenadas  de  las  cel¬ 
das.  Las  transformaciones  geométricas,  por  el  contra¬ 
rio,  conservan  la  radiometría  pero  modifican  la  geo¬ 
metría  haciendo  variar  la  posición  de  los  puntos  de 
la  imagen. 


Estadísticos  básicos  de  la  imagen 


Para  el  tratamiento  digital  de  la  imagen  va  a  ser  re¬ 
currente  echar  mano  de  algunos  estadísticos  básicos 
de  la  imagen.  Para  ello  se  definen  los  siguientes  con¬ 
ceptos: 

Histograma 

El  histograma  de  la  imagen  es  la  representación  del 
número  de  celdas  que  tienen  cada  uno  de  los  k  nive¬ 
les  digitales  posibles.  Es  el  equivalente  estadístico  del 
histograma  de  frecuencias.  Por  ello  también  es  deno¬ 
minado  a  veces  función  de  densidad  (fdd)  o  p{z). 

Puede  asimilarse  a  una  función  probabilística  si  se 
considera  el  universo  muestral  el  conjunto  de  todas 
las  celdas  de  la  imagen,  cada  una  de  ellas  con  un 
cierto  ND.  La  frecuencia  relativa  vendría  a  ser  en  este 
caso  la  probabilidad  de  encontrar  en  toda  la  imagen 
celdas  con  un  determinado  ND. 

Histograma  acumulativo 

Es  el  histograma  del  número  de  celdas  que  en  la  ima¬ 
gen  tienen  un  determinado  nivel  digital  o  inferior  a 


Histograma 


Núm.  celdas  Suma  acumulada 


Representación  del  histograma  de  una  imagen  y  de  su  histograma  acumulativo. 
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celdas  acumuladas  (millones) 
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éste.  Coincide  con  el  concepto  estadistico  de  función 
de  distribución  (ffD)  o  P{z). 

Media 


Varianza 

Es  una  medida  de  dispersión  de  los  ND  en  torno  al 
nivel  promedio.  Encierra  una  forma  de  considerar  el 
contraste  y  la  información  que  porta  una  imagen. 


La  media  de  la  imagen  es  el  valor  promedio  de  los 
ND  del  conjunto  de  celdas. 


ND 

N 


--‘Y. 

ij 

Valores  extremos 


(nd-nd) 

N 


donde  ND  es  el  nivel  digital  promedio,  i,j  las  colum¬ 
nas  y  líneas  de  la  imagen  y  fV  el  número  de  celdas, 
esto  es  i-j. 


Será  interesante  con  diversos  propósitos  considerar 
los  niveles  digital  mínimo  NDmin  y  máximo  NDmax  de 
una  imagen,  es  decir,  el  nivel  digital  de  la  celda  más 
oscura  y  el  de  la  más  clara  de  la  imagen,  respectiva¬ 
mente,  que  en  términos  generales  no  ha  de  ser  nece¬ 
sariamente  negro  (0)  o  blanco  (^-1). 


Transformaciones  radiométricas 


Las  transformaciones  radiométricas  son  utilizadas 
para  alterar  el  brillo  de  la  imagen,  mediante  operacio¬ 
nes  de  realce  o  reducción  del  contraste.  También  pue¬ 
den  emplearse  para  minimizar,  que  no  corregir  por 
completo,  determinadas  anomalías  sufridas  durante 
la  adquisición  o  grabación  de  la  imagen. 

Se  entiende  por  rango  o  margen  dinámico  la  ampli¬ 
tud  completa  de  niveles  digitales  que  una  imagen 
puede  exhibir,  según  sea  su  resolución  radiométrica. 

Para  un  rango  dinámico  k-. 

ND  =  f{x,y)-,  Q<ND<k-\ 


Siendo  0  el  negro  y  ^-1  el  blanco  de  la  celda. 

En  el  análisis  digital  puede  ser  interesante  operar  con 
la  imagen  normalizada,  esto  es,  acotada  entre  0  y  1: 

z  =  g{x,y) 

de  modo  que: 

A^Z)  =  RED(ytz'-l) 

Siendo  RED  la  función  de  redondeo.  Por  ejemplo,  si 
z  =  0,5: 


Teledetección:  principios  y  aplicaciones 


08-6 


Tratamiento  digital  de  imágenes 


ND  =  RED[156(),5-\]  =  \11 


Redondeo  y  truncamiento 

Entre  las  transformaciones  radiométricas  más  senci¬ 
llas  se  encuentra  el  redondeo  y  truncamiento.  Si  la 
imagen  original  posee  un  rango  dinámico  0  <  ND  < 
/c  —  1  y  la  imagen  final,  que  no  tiene  por  qué  conser¬ 
varlo,  tiene  un  rango  0  <  ND'  <k'  —  1,  entonces: 


0 

ND'  <  0 

ND'  =  < 

REDfA^D] 

>  con  < 

0<ND'<k'-l 

k'-\ 

L  J 

k'-\<ND' 

El  operador  acota  entre  0  y  ^'-1  la  imagen,  aun  en  el 
supuesto  de  que,  como  resultado  de  una  transforma¬ 
ción  previa,  los  valores  pudieran  ser  negativos  o  ex¬ 
ceder  del  nuevo  rango  dinámico,  y  además  redondea 
a  valores  enteros  cualquier  nivel  de  brillo  compren¬ 
dido  dentro  de  ese  rango  dinámico. 

En  las  transformaciones  radiométricas,  la  condición 
necesaria  para  que  no  se  produzca  una  inversión  del 
contraste  y  para  que  exista  la  transformación  inversa 
es  que  T{z)  sea  biunivoca  y  monótonamente  cre¬ 
ciente,  es  decir  que: 

r(zi)>r(z2)  si  Zi>z2 

para  lo  cual  es  necesario  que  el  rango  dinámico  de  la 
imagen  final  no  sea  inferior  al  de  la  imagen  original: 


Carlos  Pinilla 
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Transformaciones  generales  de  amplitud 

Existen  numerosas  transformaciones  radiométricas, 
cada  una  de  ellas  orientada  a  conseguir  un  determi¬ 
nado  efecto.  Considerando  el  caso  más  general  en  el 
que  la  imagen  original  y  la  transformada  tengan,  res¬ 
pectivamente,  rangos  dinámicos  k  y  k\  pueden  ci¬ 
tarse  como  unas  de  las  más  habituales  la  identidad  y 
el  negativo: 

-1  Identidad: 


f 

z  =z 

lo  que  implica  que: 


ND'  =  RED 


k'-\ 

k-\ 


■ND 


-1  Negativo,  en  el  sentido  de  negativo  fotográfico, 
esto  es,  inversión  del  contraste: 

z  =  \-z 

lo  que  implica  que: 


ND'  =  RED 


(t-'-i/i 

V 


ND^ 


Otras  transformaciones  de  amplitud  podrían  implicar 
funciones  exponenciales  o  logarítmicas. 


Correcciones  radiométricas 


Las  correcciones  radiométricas  son  operaciones  en¬ 
caminadas  a  eliminar  o  minimizar  las  alteraciones 
que  la  imagen  original  haya  podido  sufrir  en  sus  va¬ 
lores  de  luminancia,  bien  a  nivel  general  o  bien  local¬ 
mente  como  consecuencia  de  algún  fallo  en  el  pro¬ 
ceso  de  captura  y/ o  grabación. 

Restauración  de  celdas  o  líneas  perdidas 

Durante  el  transcurso  de  la  adquisición  o  grabación 
de  la  imagen,  es  posible  que  por  diversas  causas  al¬ 
gunas  celdas  aisladas  que  se  queden  sin  un  ND  asig¬ 
nado,  presentando  valores  claramente  erróneos  y 
apareciendo  en  la  imagen  como  completamente  ne¬ 
gras  {ND  =  0)  o  completamente  blancas  {ND  =  k-\). 
Existen  técnicas  para  la  restauración  de  estas  celdas 
que  operan  en  el  dominio  de  la  frecuencia  y  que  se¬ 
rán  abordadas  en  el  capítulo  correspondiente.  En  el 
dominio  espacial  es  muy  fácil  buscar  un  remedio  a 
esta  situación  asignando  a  la  celda  perdida  en  la  po¬ 
sición  de  columna  y  fila  {i,j)  un  valor  igual  al  prome¬ 
dio  de  las  cuatro  vecinas  más  próximas: 


ND'..  =  RED 


O  al  promedio  de  las  ocho  celdas  circundantes: 
ND'  =REdÍ 

g 


{ND,,-,  +ND„,^_,  +ND„,^+ND„,„  + 

+ m.,„ + + m-,., + m-,,-, ) 


En  otras  ocasiones,  en  lugar  de  celdas  aisladas,  se 
pierden  o  se  alteran  líneas  completas  de  imagen.  Esta 


Correcciones  radiométricas 
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circunstancia  se  hace  patente  cuando  el  valor  prome¬ 
dio  de  los  ND  de  las  celdas  de  toda  la  línea  alterada 
es  significativamente  mayor  o  menor  que  el  promedio 
del  resto  de  ND  de  la  imagen.  Una  primera  forma  de 
asignar  un  valor  a  la  la  línea  perdida  es  tomar  para 
cada  celda  el  ND  de  la  celda  vecina  de  la  línea  ante¬ 
rior: 


ND',.=ND,._„ 

de  la  siguiente: 

ND',.=ND,.,„ 

o,  mejor  aún,  del  promedio  de  ambas: 

WD',=REDri(A'D,,,_,+JVD,,.,) 


En  lugar  de  las  líneas,  pueden  usarse  de  igual  modo 
las  columnas. 

Esta  operación  no  recupera  los  niveles  digitales  defi¬ 
cientemente  adquiridos,  sino  que  solamente  mini¬ 
miza  el  efecto  óptico  o  mejora  en  cierta  medida  el 
tratamiento  digital  posterior.  Además,  es  condición 
imprescindible  para  que  su  aplicación  esté  justificada 
desde  el  punto  de  vista  matemático  que  exista  alta 
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autocorrelación'  entre  los  ND  de  las  celdas  vecinas. 

Hay  veces  que  el  promedio  entre  las  líneas  anterior  y 
posterior  no  asegura  un  resultado  óptimo,  y  esto 
puede  asegurarse  a  la  vista  de  lo  que  sucede  en  una 
banda  de  la  imagen  que  esté  altamente  correlacio¬ 
nada  con  la  que  presenta  el  defecto  {K).  En  esta 
banda  de  referencia  {R)  podría  suceder  que  el  prome¬ 
dio  de  las  líneas  próximas  difiera  ligeramente  de  los 
valores  promedio  de  la  línea  central  con  la  que  se 
compara.  Se  puede  trasladar  esta  diferencia  a  la 
banda  defectuosa  y  realizar  una  corrección  más  de¬ 
purada: 


NDf:^.  =  RED 


j+i 


-I- 


de  referencia  para  homologarlos  a  los  de  la  banda 
defectuosa. 

Bandeado 

Para  optimizar  el  rendimiento  de  los  sensores  basa¬ 
dos  en  sistemas  de  escáner,  con  cada  oscilación  del 
espejo  no  se  rastrea  una  sola  línea  sino  un  bloque  de 
varias.  Para  ello,  en  el  plano  focal  se  disponen  tantos 
detectores  como  líneas  constituyan  dicho  bloque.  Si 
el  detector  n-simo  se  descalibra,  todas  las  líneas  ba¬ 
rridas  por  él  quedarán  registradas  con  niveles  anó¬ 
malos  de  brillo,  bien  por  defecto  o  bien  por  exceso. 
Como  el  fallo  es  permanente,  la  imagen  exhibirá  lí¬ 
neas  anómalas  en  intervalos  regulares  de  n  líneas 
Esta  distorsión  se  conoce  como  bandeado  de  la  ima¬ 
gen. 


-t- 


ND, 


ij 


donde  los  superíndices  {K  y  {R  hacen  referencia  a  las 
bandas  a  corregir  y  la  de  referencia,  respectivamente 
y  s  es  la  desviación  típica  en  cada  una  de  las  bandas. 
El  último  término  de  la  derecha  representa  el  prome¬ 
dio  de  los  ND  de  las  líneas  en  la  banda  R  y  todo  lo 
que  está  dentro  del  paréntesis  es  la  diferencia  entre 
el  valor  real  de  la  línea  central  y  el  promedio  de  las 
circundantes.  El  primer  término  dentro  del  corchete 
es  el  promedio  de  ND  en  la  banda  Ky  ú  cociente  de 
las  desviaciones  típicas  es  el  factor  de  expansión  que 
debe  aplicarse  a  los  resultados  tomados  de  la  banda 


Es  posible  detectar  si  hay  bandeado  construyendo  y 
comparando  entre  sí  los  histogramas  parciales  de  la 
imagen  que  corresponden  al  conjunto  de  líneas  regis¬ 
tradas  con  cada  detector.  En  condiciones  normales 
todos  los  histogramas  tendrán  estadísticos  similares, 
particularmente  los  valores  de  nivel  digital  medio  des¬ 
viación  típica.  Pero  si  un  detector  falla,  es  probable 
que  todas  las  líneas  por  él  registradas  presenten  ni¬ 
veles  digitales  inferiores,  por  ejemplo,  y  una  menor 
varianza  o  desviación  típica.  La  corrección  del  ban¬ 
deado  consiste  en  conseguir  para  el  conjunto  de  lí¬ 
neas  barridas  por  ese  detector  los  mismos  estadísti¬ 
cos  que  para  el  resto  de  la  imagen.  Si  k  es  el  conjunto 
de  las  líneas  defectuosas: 


*  La  autocorrelación  es  una  función  estadística  definida  como  la  correlación 
cruzada  de  una  señal,  en  este  caso  una  imagen,  consigo  misma.  La  autoco¬ 
rrelación  mide  el  grado  de  similitud  entre  una  función  y  una  versión  de 
ésta  desplazada.  En  el  caso  de  la  imagen,  la  autocorrelación  de  los  ND  viene 
a  significar  que  el  valor  de  una  celda  depende  estrechamente  de  los  valores 
de  las  celdas  vecinas,  al  menos  dentro  de  un  contexto  determinado.  Los 


modelos  digitales  del  terreno  (MDT)  son  funciones  que  tienen  gran  autoco¬ 
rrelación  dentro  del  contexto  de  una  ladera,  limitada  por  las  líneas  de  ro¬ 
tura.  Para  la  mayoría  de  los  terrenos  de  topografía  suave,  la  cota  en  un 
punto  depende  mucho  de  la  cota  en  puntos  cercanos.  Con  la  imagen  de 
reflectancia,  ocurre  algo  similar:  es  frecuente  que  los  valores  de  luminancia 
no  cambien  bruscamente  de  una  celda  a  la  vecina. 


Imagen  Histogramas  de  cada  detector 


Esquema  de  proceso  de  la  corrección  del  bandeado 
de  una  imagen  adquirida  por  un  sensor  de  barrido 
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plano  de  barrido 
línea  de  barrido 


traza 


troyectorb  de  vuelo 


plano  de  iluminación 
celdilla  en  captura 


NDlj  =  RED 


(^ND.j  +  NDk)  +  ND 


donde  s  es  la  desviación  típica  de  los  ND  de  toda  la 
imagen,  a  excepción  del  subconjunto  defectuoso, 

cuya  desviación  típica  es  Sk  y  ND  y  NDk  son  los  ni¬ 
veles  digitales  promedio  de  la  imagen,  a  excepción 
del  conjunto  de  líneas  defectuosas  y  el  promedio  de 
estas  últimas,  respectivamente. 

Si  se  escribe: 


AíD;j=RED[a./VD,j+í..] 


y  además  las  celdas  son  vistas  con  diferentes  grados 
de  inclinación. 

En  definitiva,  interesa  estudiar  el  triedro  que  forman 
el  vector  característico  perpendicular  a  la  celda,  la  vi¬ 
sual  instantánea  desde  el  sensor  y  el  vector  de  orien¬ 
tación  solar  para  corregir  la  radiancia  simulando  una 
superficie  plana  horizontal  y  una  visual  ortogonal  al 
terreno.  La  radiancia  estimada  como  real  puede  ser 
deducida  a  partir  de  la  que  ha  registrado  la  imagen 
del  siguiente  modo: 


Vectores  y  ángulos  de  interés  en  la  adquisición  de  la 
imagen:  n  =  vector  característico;  v  =  visual  instantá¬ 
nea.;  s  =  orientación  solar;  a  =  inclinación  del  terreno; 
K  =  orientación  del  terreno;  cp  =  acimut  solar;  ^  =  eleva¬ 
ción  solar;  0  =  visión  instantánea;  s  =  inclinación  de  la 
órbita;  i  =  incidencia  solar  sobre  el  terreno;  e  =  inclina¬ 
ción  de  la  observación. 


puede  identificarse  los  coeficientes  a  y  b  con  la  ga¬ 
nancia  y  sesgo,  respectivamente: 

a,=—;  b,=ÑD-—ÑDk 

Normalización  topográfica 

Con  el  nombre  de  normalización  topográfica  se  cono¬ 
cen  una  serie  de  procedimientos  encaminados  a  co¬ 
rregir  el  efecto  que  la  topografía  de  la  superficie  in¬ 
duce  en  su  reflectancia  aparente.  Es  sabido  La  radian¬ 
cia  mantiene  una  relación  lineal  con  el  nivel  digital 
de  la  celda  y  que  lo  que  ciertamente  caracteriza  a  la 
celda  es  su  reflectancia,  la  cual  también  mantiene  la 
linealidad  con  el  nivel  digital  cuando  la  superficie  es 
horizontal.  Sin  embargo,  en  condiciones  de  orografía 
muy  acusada,  la  radiancia  direccional  no  representa 
suficientemente  bien  a  la  reflectancia,  puesto  que 
aquélla  depende  de  la  irradiancia  y  del  ángulo  de  vi¬ 
sión  del  sensor.  Hay  diferencias  de  iluminación  entre 
unas  laderas  y  otras  dependiendo  de  su  orientación. 


en  cuya  expresión,  es  la  radiancia  estimada  como 
real  procedente  de  la  superficie  plana,  es  la  ra¬ 
diancia  registrada  en  la  imagen,  T(p,q)  es  el  factor  de 
normalización  topográfica  para  una  celda  de  pendien¬ 
tes  py  q,  L{0,0)  es  la  radiancia  de  la  superficie  hori¬ 
zontal  y  L{p,q)  es  la  radiancia  de  la  superficie  incli¬ 
nada. 

Existen  diferentes  procedimientos  aplicara  la  norma¬ 
lización  topográfica,  algunos  de  los  cuales  se  mues¬ 
tran  a  continuación: 

— 1  Normalización  por  el  coseno  del  ángulo  de  inci¬ 
dencia  solar: 

.  COS0, 

T  = - ^ 

cosí 

siendo  0s  el  ángulo  cenital  solar  e  i  el  ángulo  de 
incidencia  solar  sobre  la  celda.  Este  método  ado¬ 
lece  del  defecto  de  sobrecorregir  las  imágenes 
cuando  eos/  es  muy  pequeño. 
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Normalización  de  C-Huang  Wei  para  cubiertas 
lambertianas: 


L  =  L-L,,„ 


^cose^-cos¿^-„^ 

.  .  ^mín 

COSí-COSí„,„  j 


donde  Lmin  es  la  radiancia  mínima  antes  de  la  co¬ 
rrección. 


empírica. 


Corrección  por  superficies  de  Minnaert: 

í  ^ 


L  =  L 


cose 


vcos^  /-cos^  e  j 


— 1  Corrección  basada  en  la  geometría  sol-cubierta- 
sensor: 

-  ,  f  COS0  +cosS-i-C^ 

L  =  L-  - - 

eos  i  +  C  J 

donde  S  es  la  pendiente  local  y  C  es  una  constante 


Esta  corrección  introduce  los  ángulos  de  inciden¬ 
cia  solar  i  y  el  de  visión  de  la  celda  desde  el  sensor 
e,  afectados  por  un  coeficiente  k  de  Minnaert,  que 
tiene  en  cuenta  el  carácter  lambertiano  de  la  su¬ 
perficie  y  que  oscila  entre  0,  para  superficies  abso¬ 
lutamente  lambertianas  y  1,  para  superficies  espe¬ 
culares. 


Correcciones  atmosféricas 


Los  aerosoles  atmosféricos  causan  un  efecto  de  dis¬ 
persión  que,  si  bien  pueden  absorber  radiación  en  de¬ 
terminados  intervalos  de  longitudes  de  onda,  común¬ 
mente  incrementan  la  irradiancia  sobre  el  sensor  en 
el  espectro  solar.  Este  efecto  de  la  dispersión  no  es 
constante  ni  en  el  tiempo,  ni  en  el  espacio.  Por  esta 
razón,  las  correcciones  generales  sólo  pueden  ser 
consideradas  como  aproximaciones,  en  tanto  que 
para  una  corrección  rigurosa  deberán  realizarse  me¬ 
didas  en  tierra  simultáneas  a  la  adquisición  de  la  ima¬ 
gen  y  aplicar  correcciones  acordes  con  el  efecto  su¬ 
frido,  por  lo  general  diferentes  para  cada  celda  de  la 
imagen. 

La  corrección  del  efecto  atmosférico  es  imprescindi¬ 
ble  cuando  se  desee  deducir  cualquier  magnitud  fí¬ 
sica  a  partir  de  los  ND  de  la  imagen,  cuando  hayan 
de  efectuarse  combinaciones  lineales  entre  bandas, 
por  ser  la  dispersión  un  fenómeno  dependiente  de  la 


longitud  de  onda  y  cuando  se  realicen  estudios  mul- 
titemporales,  en  cada  una  de  cuyas  fechas,  la  reper¬ 
cusión  déla  atmósfera  puede  haber  sido  diferente. 

Componentes  de  la  radiancia  ascendente 

Si  se  considera  un  sensor  con  un  campo  de  visión 
instantáneo  suficientemente  estrecho  en  relación  con 
la  variabilidad  de  la  superficie  explorada,  la  radiancia 
capturada  por  el  instrumento  puede  descomponerse 
en  una  suma  de  radiancias  de  diversa  procedencia. 

En  primer  lugar,  como  ya  se  dijo  en  el  capítulo  4,  se 
encuentra  la  señal  Ls  atenuada  producida  por  la  ra¬ 
diancia  de  la  luz  que  atraviesa  la  atmósfera  en  sen¬ 
tido  descendente,  que  se  refleja  en  la  superficie  y  que 
es  transmitida  hacia  el  sensor  de  nuevo  a  través  de  la 
atmósfera.  De  entre  todas  las  demás,  ésta  es  la  señal 
útil  que  sirve  para  identificar  la  superficie  y  obtener 


Relación  entre  la  situación  del  sensor 
y  el  ángulo  de  visión  de  la  celda. 


Carlos  Pinilla 


de  ella  sus  características  físicas. 


Desgraciadamente,  no  es  la  única  señal  que  incide  en 
el  plano  focal  del  sensor.  También  penetra  radiancia 
directa  Lo  que  procede  del  sol  y  que,  previa  disper¬ 
sión  en  el  seno  de  la  atmósfera  y  sin  haber  alcanzado 
la  superficie,  termina  incrementando  la  señal. 

Y  todavía  hay  dos  fuentes  más  de  radiancia  indirecta. 
De  una  parte,  la  luz  Ld\  que  ha  sido  reflejada  en  la 
dirección  del  sensor  por  una  superficie  que  no  ha 
sido  iluminada  directamente  por  el  sol,  sino  a  partir 
de  la  luz  difusa  atmosférica,  siendo  el  caso  más  ex¬ 
tremo  es  el  de  las  superficies  en  sombra,  y  finalmente, 
puede  penetrar  en  el  instrumento  radiancia  Ld2  refle¬ 
jada  por  una  superficie  que  ha  sido  iluminada  direc¬ 
tamente  por  el  sol  o  por  la  propia  atmósfera,  pero 
cuya  luz  es  dispersada  al  menos  una  vez  dentro  de  la 
atmósfera  antes  de  ser  captada  por  el  sensor.  Dicho 
de  otra  forma,  se  trata  de  celdas  iluminadas  directa  o 
indirectamente  que  no  están  bajo  el  IFOV  del  sensor. 

En  definitiva,  se  tiene  que: 

L  =  Lq+L^  +  L¿j  -i-  L¿2 

y  las  correcciones  atmosféricas  que  se  implementen 
deberán  ir  encaminadas  a  estimar  de  la  forma  más 
aproximada  el  valor  de  Ls  a  partir  del  valor  medido 
L. 

Corrección  por  los  valores  mínimos 

Una  de  las  correcciones  más  rápidas  y  extendidas  es 
la  denominada  corrección  del  histograma  por  sus 
valores  mínimos.  Este  procedimiento  asume  que  el 
ND  mínimo  en  los  canales  visibles,  especialmente  en 
el  rojo,  debe  corresponder  a  las  zonas  de  agua  limpia 
y  profunda.  A  la  vista  de  la  signatura  espectral  de  esta 
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superficie,  en  ausencia  de  efecto  atmosférico  la  reflec- 
tancia  debería  ser  próxima  a  cero.  En  consecuencia, 
los  valores  mínimos  distintos  de  cero  de  los  histogra- 
mas  de  las  bandas  afectadas  han  de  estar  causados 
por  la  radiación  difusa  atmosférica.  El  método  de  co¬ 
rrección  consiste  en  restar  de  todos  los  ND  de  cada 
banda  su  ND  mínimo.  Para  la  banda  k,  por  ejemplo: 


ND'  ,,,=ND.  ,,,-ND  . 

i,j{k)  i,j{k)  mm(k) 


Es  importante  tener  en  cuenta  que  la  consideración 
anterior  se  refiere  a  superficies  de  agua  con  las  ca¬ 
racterísticas  citadas,  las  cuales  no  siempre  están  dis¬ 
ponibles  en  la  imagen.  Es  relativamente  frecuente  en¬ 
contrar  cuerpos  de  agua  más  o  menos  profunda, 
como  la  retenida  en  embalses  o  lagos,  que,  sin  em¬ 
bargo,  no  puede  considerarse  limpia  a  los  efectos 
que  interesan  aquí.  A  menudo  son  aguas  eutrofizadas 
y  ricas  en  clorofila,  cuya  radiancia  en  el  visible  no  es 
tan  baja  como  la  del  agua  marina.  En  el  caso  de  los 
ríos  o  cabeceras  de  embalses  hay  que  advertir  tam¬ 
bién  la  presencia  de  fondos  poco  profundos  que  con¬ 
tribuyen  en  la  reflexión  de  la  luz  solar.  En  estos  casos 
hay  que  ser  prudente  a  la  hora  de  restar  los  valores 
mínimos  para  no  incurrir  en  una  sobrecorrección.  Si 
no  se  dispone  de  más  datos,  podría  considerarse  la 
posibilidad  de  corregir  sustrayendo  tan  solo  de 
los  niveles  digitales  mínimos. 


Correcciones  sobre  superficies  de  tierra 

En  numerosas  ocasiones  puede  no  encontrarse  en  la 
imagen  una  sola  superficie  de  agua  apta  para  utilizar 
sus  ND  para  la  corrección  atmosférica.  En  estos  casos 
debe  recurrirse  a  zonas  de  baja  radiancia,  como  por 
ejemplo  ciertas  manchas  de  bosque  denso  y  som¬ 
breado,  delimitarlas  sobre  la  escena  y  extraer  los  ND 
en  las  distintas  bandas.  A  partir  de  estos  valores  se 
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calcula  la  recta  de  regresión  entre  cada  banda  y  al¬ 
guna  otra  de  referencia  en  la  que  la  dispersión  atmos¬ 
férica  tenga  menor  repercusión,  por  ejemplo  el  rojo. 

La  ordenada  en  el  origen  de  la  recta  de  regresión  in¬ 
dica  el  valor  a  sustraer  de  todos  los  ND  de  esa  banda. 

La  vegetación  densa  y  pobremente  iluminada,  como 
es  la  de  los  bosques  situados  en  laderas  orientadas  al 
norte  en  el  hemisferio  boreal  y  en  ladera  sur  en  el 
austral,  tiene  baja  reflectancia  en  el  visible.  Esta  cua¬ 
lidad  es  ampliamente  validable  para  bosques  en  cual¬ 
quier  latitud  y  época,  salvo,  por  motivos  evidentes, 
los  caducifolios  en  invierno,  y  se  encuentra  recogida 
en  la  literatura  especializada.  Para  ángulos  cenitales 
solares  próximos  a  40”  la  reflectancia  de  las  praderas 
oscila  entre  el  2,6  %  en  el  azul  y  el  3,9  %  en  el  rojo, 
en  tanto  que  para  los  bosques  caducifolios  está  entre 
el  1  y  el  2,5  %  en  ambos  canales.  Cualquier  cálculo 
de  reflectancia  sobre  superficies  semejantes  que  se 
aleje  de  estas  cifras  deberá  ser  entendido  como  pre¬ 
sencia  de  efecto  atmosférico. 

Para  la  aplicación  del  método  solamente  es  necesaria 
una  imagen  que  contenga  alguna  zona  de  vegetación 
densa  y  sombreada. 

Corrección  atmosférica  en  el  análisis 
muititemporai 

Los  estudios  multitemporales  exigen  comparar  imá¬ 
genes  con  similar  influencia  atmosférica,  dado  que 
cualquier  diferencia  en  la  radiometría  de  las  imáge¬ 
nes  imputable  a  efecto  atmosférico  podría  ser  inter¬ 
pretada  como  un  cambio,  aun  cuando  en  realidad  la 
superficie  pueda  haber  permanecido  inmutable. 

La  corrección  puede  hacerse  de  modo  relativo  deter¬ 
minando  la  reflectancia  de  una  serie  de  cubiertas  de 
reflectancia  invariable  a  lo  largo  del  año  y  calculando 


los  coeficientes  lineales  de  ganancia  y  sesgo  de  esas 
cubiertas  en  cualquier  fecha  con  respecto  a  la  escena 
de  referencia.  Para  normalizar  radiométricamente  la 
reflectancia  a  los  valores  de  la  fecha  de  referencia  se 
aplica  la  transformación  inversa  con  esos  mismos 
coeficientes  para  todas  las  celdas  de  la  imagen. 

En  el  capítulo  dedicado  la  teledetección  multitempo- 
ral  se  dan  más  detalles  de  los  procedimientos  a  em¬ 
plear. 

Corrección  atmosférica  rigurosa 

La  corrección  atmosférica  rigurosa  requiere  tener  un 
conocimiento  preciso  de  los  parámetros  atmosféricos 
sobre  la  escena  en  el  momento  de  la  captura  de  la 
imagen.  Los  más  importantes  son  el  coeficiente  de 
extinción,  el  espesor  óptico,  el  albedo  por  dispersión 
simple  y  la  función  de  fase  de  dispersión,  y  pueden 
obtenerse  a  partir  de  distintas  fuentes: 

-1  Climatología.  Se  pueden  hacer  estimaciones  para 
un  determinado  lugar  y  fecha  a  partir  de  bancos 
de  datos  tabulados  de  las  características  atmosfé¬ 
ricas  para  cada  estación  del  año. 

-1  Medidas  de  campo.  Se  miden  los  parámetros  at¬ 
mosféricos  mediante  observaciones  fotométricas 
direccionales  del  sol  y  por  recolección  de  partícu¬ 
las  en  suspensión  en  las  altas  capas  de  la  atmós¬ 
fera  con  aviones  o  globos. 

-1  Medidas  sobre  la  imagen.  El  albedo  por  dispersión 
simple  y  la  función  de  fase  de  dispersión,  pueden 
medirse  sobre  imágenes  con  líneas  de  fuerte  con¬ 
traste,  como  por  ejemplo  son  los  contactos  tierra- 
agua. 

Existen  modelos  de  corrección  atmosférica  muy  pre¬ 
cisos  basados  en  algoritmos  que  comparan  la  reflec¬ 
tancia  de  zonas  de  alto  contraste. 


Banda K 


Banda R 


Corrección  por  regresión  sobre  otra  banda 
poco  influida  por  el  efecto  atmosférico. 


Carlos  Pinilla 
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Realce  de  imágenes 

Se  conocen  como  operaciones  de  realce  de  imágenes 
las  diversas  técnicas  orientadas  a  mejorar  o  acentuar 
la  calidad  de  las  imágenes  con  el  objetivo  de  facilitar 
el  análisis  visual,  sobre  papel  o  monitor,  o  para  per¬ 
feccionar  la  producción  cartográfica. 

Se  puede  hablar  de  dos  tipos  de  realce,  el  radiomé- 
trico,  basado  en  tareas  de  expansión  del  contraste,  y 
el  geométrico,  sustentado  por  el  filtrado  digital. 


Operaciones  con  el  contraste 

El  contraste  es  un  indicativo  de  la  amplitud  del  inter¬ 
valo  de  ND  presentes  en  una  imagen,  de  entre  la 
gama  completa  de  valores  que  podría  ofrecer. 

Existen  varias  formas  de  medir  el  contraste  en  una 
determinada  banda  de  la  imagen; 


C2=ND^^-ND^^ 


^3  ~  ^ND 


La  primera  de  ellas,  el  cociente  entre  el  ND  máximo 
y  el  mínimo  da  una  idea  relativa  del  rango  dinámico 
de  la  imagen.  La  diferencia  proporciona  la  misma  in¬ 
formación,  pero  desde  un  punto  de  vista  absoluto.  Fi¬ 
nalmente,  la  varianza  del  histograma  de  la  imagen  in¬ 
forma  de  la  dispersión  de  los  valores  a  partir  del  valor 
medio.  Entendido  de  este  modo,  la  varianza  (o  en  su 
caso,  la  desviación  típica)  no  sólo  es  una  medida  del 


contraste  sino  también  de  la  cantidad  de  información 
contenida  en  la  imagen. 

Para  una  determinada  resolución  radiométrica,  la 
imagen  puede  ocupar  todo  el  rango  dinámico,  o  no. 
En  el  primer  caso,  NDmin=0  y  NDmax=k-\,  siendo  k  el 
número  de  niveles  digitales  representables.  En  una 
gama  acromática,  ambos  valores  serán  representados 
con  un  color  negro  y  blanco,  respectivamente. 

Sin  embargo,  por  diversas  razones,  la  imagen  puede 
no  abarcar  todo  el  rango  dinámico  del  sistema.  Tén¬ 
gase  en  cuenta  que,  al  diseñar  el  sensor,  debe  estar 
preparado  para  adquirir  imágenes  en  lugares  muy  de 
la  superficie  terrestre  dispares,  desde  las  de  muy  alta 
reflectancia,  como  los  casquetes  polares  en  el  visible, 
hasta  otras  de  reflectancia  prácticamente  nula,  como 
los  océanos  en  el  infrarrojo  próximo  o  la  vegetación 
umbrosa  en  el  rojo  espectral.  Si  en  una  escena  parti¬ 
cular  no  existen  esas  superficies  extremas,  la  gama 
de  ND  exhibida  no  ocupará  todo  el  rango  dinámico, 
sino  sólo  un  conjunto  de  valores  intermedios.  En  con¬ 
secuencia,  la  celda  más  oscura  de  la  imagen  mono¬ 
croma  no  será  negra,  sino  un  gris  más  o  menos  os¬ 
curo,  y  la  más  brillante  tampoco  será  blanca  sino  un 
valor  de  gris  más  o  menos  claro. 

Las  operaciones  de  ajuste  del  contraste  pueden  ser 
realizadas  sobre  la  imagen  directamente,  generando 
una  nueva  imagen  de  salida  con  los  ND  alterados,  o 
bien  pueden  ser  realizadas  sobre  la  pantalla,  simple¬ 
mente  alterando  el  brillo  que  corresponde  a  cada  ND, 
pero  sin  modificar  el  archivo  original  de  la  imagen  y 
sin  generar  una  nueva. 
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Se  define  el  nivel  visual  {NV)  como  el  brillo  atribuido 
a  una  celda  de  imagen  en  el  sistema  de  visualización 
o  de  impresión.  Siempre  que  no  haya  mediado  una 
alteración  del  contraste,  ND=NV. 

La  relación  entre  ND  y  NV  puede  representarse  me¬ 
diante  una  gráfica  en  la  que  el  ND  se  sitúa  en  absci¬ 
sas  y  el  NV  en  ordenadas.  El  caso  anterior  de  ausen¬ 
cia  de  alteración  viene  representado  por  una  recta  de 
pendiente  unidad.  Cualquier  operación  que  lleve  apa¬ 
rejada  una  reducción  del  contraste  implicará  una  re¬ 
lación  NV-ND  lineal  de  pendiente  <  1.  Por  el  contra¬ 
rio,  la  expansión  del  contraste  vendrá  establecida  me¬ 
diante  una  función,  en  este  caso  lineal,  de  pen¬ 
diente  >  1. 

En  general  no  es  interesante  reducir  el  contraste  de 
una  imagen,  dado  que  esta  operación  siempre  lleva 
consigo  una  pérdida  de  información.  Solamente  es 
imprescindible  realizarlo  cuando  el  sistema  de  trata¬ 
miento  o  de  visualización  tiene  un  rango  dinámico 
inferior  a  de  la  propia  imagen.  En  este  caso  hay  que 
reunir  bajo  un  mismo  nivel  digital  nuevo  rangos  de 
valores  originales  que  implicaban  más  de  un  nivel  di¬ 
gital.  Obviamente,  todas  las  celdas  de  brillo  dentro  de 
ese  rango  que  en  la  imagen  original  eran  distingui¬ 
bles  pasarán  ahora  a  ser  iguales. 

Las  operaciones  más  frecuentes  de  alteración  del  con¬ 
traste,  pues,  son  las  de  expansión.  Se  aplican  estas 
transformaciones  para  adecuar  la  resolución  radio- 
métrica  de  la  imagen  a  la  capacidad  de  visualización 
del  sistema,  representada  por  el  rango  dinámico  de  la 
tarjeta  gráfica  o  del  número  de  niveles  de  gris  que 
pueda  proporcionar  el  periférico  de  salida.  Así,  me¬ 
diante  la  expansión  del  contraste  de  la  imagen  se 
hace  corresponder  un  nuevo  rango  de  ND  mayor,  o 
de  NV  si  sólo  se  pretende  alterar  la  muestra  en  pan¬ 
talla,  que  le  rango  de  ND  presentes  en  la  imagen  ori¬ 
ginal.  Dado  que  el  número  de  ND  diferentes  en  ésta 


es  inferior  al  margen  dinámico  posible,  cuando  se 
haya  expandido  el  contraste  y  se  hayan  distribuido 
los  nuevos  ND  en  dicho  margen,  podrán  quedar  ND 
sin  asignar  a  las  celdas,  es  decir,  el  histograma  ten¬ 
derá  a  presentar  huecos. 


Huecos  en  el  histograma  tras  la  expansión 
del  contraste. 


Expansión  lineal 

El  más  general  de  los  casos  de  expansión  de  con¬ 
traste  es  el  lineal.  Mediante  este  tipo  de  expansión, 
en  una  gama  de  brillo  z  entre  0  y  1  se  ocupará  todo 
el  margen  dinámico  de  la  imagen  de  salida,  asignán¬ 
doles  el  valor  0  a  las  celdas  cuyo  z  sea  inferior  o  igual 
a  un  determinado  valor  umbral  mínimo  en  la  imagen 
original  {m)  y  el  valor  1  a  las  celdas  con  z  iguales  o 
superiores  a  un  determinado  valor  umbral  máximo 


Visualización  de  una  imagen  con  (a)  con¬ 
traste  original,  (b)  reducido  y  (c)  expandido. 
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frecuencia  1 


frecuencia 


0  1 

z 

Expansión  lineal  sin  saturación. 


Z 

Expansión  lineal  con  saturación. 


{M).  Entre  ambos  extremos,  se  distribuirán  lineal¬ 
mente  los  nuevos  valores  de  brillo  expandidos  {z'): 

0  z^m 

,  z-m  _  ^  ^  w 

z  =■{ -  \/m<z^M 

M  —m 

1  Vz>M 


Máxima  expansión  iineai  sin  saturación 

Con  esta  expansión  es  estira  el  histograma  de  la  ima¬ 
gen  original  hasta  hacer  cero  su  ND  mínimo  y  k’-l  su 
ND  máximo,  consiguiendo  que  la  gama  de  grises  pri¬ 
mitiva,  que  no  ocupaba  todo  el  margen  dinámico  pase 
ahora  a  discurrir  entre  el  negro  y  el  blanco.  Para  ello, 
en  la  expansión  genérica  indicada  anteriormente  se 
fijan  los  valores  m  y  71/  de  esta  forma: 


O,  si  se  prefiere: 


0 


ND'  =  \ 


ND-m 
M  -m 


ik'-l) 


k'-l 


ND<m-{k-\) 
y  m-{k-\)<ND<M  -{k-l) 
y  ND>M -{k-l) 


La  expresión  de  ND '  para  los  valores  de  ND  com¬ 
prendidos  entre  my  M  puede  descomponerse  en  dos 
sumandos,  de  modo  que  el  nivel  digital  resultante 
queda  expresado: 

,  k'-\  {k'-í)-m 

ND'  = - ND-^ - —  =  g-ND  +  s 

M—m  M  —m 


Donde  la  ganancia  {g)  y  el  sesgo  (s)  son  los  factores 
multiplicativo  y  aditivo  de  la  transformación  lineal: 

k’-\  {k'-\)-m 

g= - ;  .s  =  -^ - ^ — 

M-m  M  -  m 

La  forma  en  cómo  se  establezcan  los  valores  m  y  M 
da  lugar  a  los  distintos  tipos  de  expansión  lineal. 


m  =  ND 
M  =ND 


min 


ND'-- 


0 

ND-ND^ 


ND  -ND 

max  m 

k'-l 


V7VD<7VD^„ 

■{k'-l)  VND^^<ND<ND^ 
V  7VD  >  7VD„„ 


Una  forma  intuitiva  de  presentar  la  expansión  del 
contraste  consiste  en  construir  una  gráfica  en  la  que 
se  presenta  el  histograma  de  la  imagen  original  en 
abscisas.  Sobre  él  se  dibuja  la  función  de  expansión 
determinada  por  las  expresiones  anteriores  y  se  hace 
reñejar  la  situación  de  cada  ND  (o  z,  en  su  caso)  so¬ 
bre  ella  para  terminar  proyectándola  sobre  el  eje  au¬ 
xiliar  de  ordenadas,  en  el  cual  aparecerán  los  nuevos 
valores  de  ND’  asignados.  En  el  caso  de  la  expansión 
sin  saturación,  dado  que  la  recta  de  expansión  tiene 
su  pie  en  el  valor  mínimo  m  y  la  cabeza  en  el  valor 
M  máximo,  del  histograma  resultante  ocupará  el  mar¬ 
gen  dinámico  completo,  entre  0  y  1  si  se  habla  de 
valores  z’  o  entre  0  y  k’-\  si  se  habla  de  ND. 


Expansión  lineal  con  saturación 

Puede  darse  el  caso  de  que  una  sola  o  muy  pocas 
celdas  de  la  imagen  ocupen  ND  muy  oscuros,  incluso 
el  propio  valor  0,  bien  porque  realmente  existan  esos 
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niveles  de  brillo  real  en  la  imagen,  o  bien  porque  han 
podido  perderse  celdas  en  el  transcurso  de  la  adqui¬ 
sición  o  grabación.  Igual  puede  suceder  con  los  valo¬ 
res  altos  del  histograma.  Unas  pocas  celdas  pueden 
ser  blancas  completamente  o  tener  un  valor  muy  cer¬ 
cano  al  máximo.  En  cualquiera  de  los  casos  la  expan¬ 
sión  lineal  anterior  no  proporcionará  buenos  resulta¬ 
dos  porque  el  margen  dinámico  ya  está  ocupado,  aun¬ 
que  sea  por  muy  pocas  celdas.  Si  se  quiere  dotar  de 
mayor  margen  a  una  gran  proporción  de  la  imagen 
aplicando  una  expansión  lineal  no  queda  otra  opción 
que  escoger  otros  umbrales  my  M  diferentes  de  los 
ND  minimo  o  máximo.  Si  se  eligen  valores  del  histo¬ 
grama  que  dejen  a  su  izquierda  un  determinado  por¬ 
centaje  de  celdas  más  oscuras  y  otro  similar  de  celdas 
más  claras  a  su  derecha  y  se  expanden  linealmente 
los  valores  comprendidos  entre  a  mbos,  se  dice  que 
la  expansión  es  lineal  con  una  saturación  igual  a  ese 
porcentaje. 

Para  implementar  el  cálculo  de  la  forma  más  fácil  se 
recurre  a  construir  el  histograma  acumulativo  de  la 
imagen  y  transformarlo  en  porcentaje  sobre  el  total 
de  celdas  a  partir  del  número  de  celdas  acumuladas. 
De  este  modo,  seleccionando  la  ordenada  del  histo¬ 
grama  se  podrá  deducir  directamente  en  abscisas  el 
nivel  digital  que  deja  por  debajo  el  porcentaje  elegido 
de  celdas  más  oscuras.  De  un  modo  equivalente  se 
procede  para  las  celdas  brillantes. 

Llamando  v  a  la  proporción  de  celdas  que  se  desea 
saturar,  en  tanto  por  uno,  y  P  al  histograma  acumu¬ 
lativo  en  valores  de  tanto  por  uno,  siendo  uno  el  total 
de  celdas  de  la  imagen: 

m=z,  (^Zmin)  P{z)<v 

M  =  z,  p(z)>l-v 


en  donde  za  y  zb  son  los  valores  de  brillo  de  la  imagen, 
o  en  su  caso  los  niveles  digitales,  correspondientes  a 
las  proporciones  v  y  1-v  del  histograma  acumulativo. 

En  la  práctica  este  tipo  de  expansión  permite  ganar 
mayor  margen  dinámico  para  un  cierto  intervalo  de 
niveles  digitales  originales,  en  detrimento  de  los  va¬ 
lores  más  oscuros  o  más  claros,  que  ahora  quedarán 
saturados  como  negro  y  blanco  respectivamente. 

Expansión  con  saturación  basada  en  la 
distribución  normal 

Este  procedimiento  de  expansión  satura  los  valores 
de  gris  basándose  en  los  parámetros  de  la  distribu¬ 
ción  normal.  Asumiendo  que  el  histograma  se  ajusta 
a  este  tipo  de  distribución,  se  eligen  los  valores  m  y 
M  como  un  número  entero  de  desviaciones  típicas 
hacia  la  izquierda  y  hacia  la  derecha  del  valor  medio 
del  histograma,  respectivamente: 


Expansión  lineal  con  saturación  basada  en 
la  distribución  normal. 


m  =[i-nG 
M  =\ji  +  n<j 

Siendo  p  la  media  aritmética  de  los  ND,  <j  la  desvia¬ 
ción  típica  de  la  distribución  y  «  un  número  entero 
elegido  por  el  analista. 

Expansión  con  saturación  por  selección 
directa 

Indudablemente,  la  expansión  con  saturación  puede 
realizarse  seleccionando  directa  y  arbitrariamente  los 
valores  umbral  m  y  M,  de  tal  modo  que  se  consiga  asi 
destacar  el  rango  de  niveles  digitales  que  interese. 

Es  frecuente  necesitar  estudiar  zonas  de  la  imagen 


Expansión  lineal  con  saturación  por  selec¬ 
ción  directa  de  los  umbrales. 
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Z 


Expansión  lineal  ajustada  a  estadísticos 
prefijados. 


Expansión  mediante  ecualización  del  histo- 
grama. 


que  resulta  excesivamente  oscura  para  el  análisis  vi¬ 
sual,  como  por  ejemplo  las  superficies  marinas  o  la 
vegetación  en  el  visible.  En  estos  casos  es  útil  recu¬ 
rrir  a  esta  modalidad  de  expansión,  aunque  durante 
el  proceso  no  se  puedan  tomar  decisiones  sobre  el 
resto  de  la  imagen. 

Expansión  ajustada  a  estadísticos 

Hay  ocasiones  en  las  que  se  requiere  que  la  imagen 
adquiera  un  determinado  contraste  y  un  brillo  medio 
establecido.  Tal  es  el  caso  de  las  tareas  de  matching 
de  histogramas,  esto  es,  la  obtención  de  dos  imágenes 
con  igual  media  y  varianza  en  sus  ND,  de  modo  que 
operaciones  como  el  mosaico  puedan  ser  llevadas  a 
cabo  sin  que  se  perciba  la  desigualdad  inicial  de  los 
contrastes. 

Para  obtener  una  imagen  cuyo  histograma  presente 
una  media  y  una  varianza  elegida  con  anterioridad  se 
recurre  a  la  siguiente  transformación: 

z'  =  — (z-p)  +  p' 
a 

en  donde  p  y  a  son  el  valor  medio  y  la  varianza  de  la 
imagen  que  se  desea  ajustar  p’  y  a’  los  estadísticos 
de  referencia  a  conseguir. 

Ecualización  dei  histograma 

En  teoría  de  señal,  la  ecualización  es  una  operación 
que  se  realiza  sobre  una  señal  destinada  a  conseguir 
que  su  espectro  de  frecuencias  sea  uniforme,  es  decir, 
en  el  que  la  intensidad  de  todos  los  componentes  sea 
igual,  sea  cual  sea  la  frecuencia  básica  de  éstos.  Esto 
supone  alterar  los  coeficientes  de  Fourier  de  la  señal 
en  el  dominio  de  la  frecuencia,  o  lo  que  es  lo  mismo. 


potenciar  las  amplitudes  de  los  componentes  de  baja 
o  de  alta  frecuencia. 

Por  analogía,  en  el  procesamiento  digital  de  imáge¬ 
nes,  la  ecualización  del  histograma  de  una  imagen 
tendría  como  objetivo  obtener  una  nueva  imagen 
cuyo  histograma  tuviera  una  distribución  uniforme, 
es  decir,  que  existiera  el  mismo  número  de  celdas 
para  cada  nivel  digital.  Sin  embargo,  en  la  práctica 
esto  no  es  posible,  dado  que  la  función  de  niveles  de 
gris,  el  histograma,  es  una  función  discreta,  a  diferen¬ 
cia  de  lo  que  sucede  en  sonido. 

En  su  lugar,  y  con  el  objeto  de  que  el  balance  de  ni¬ 
veles  de  gris  sea  más  equilibrado  lo  que  se  hace  en 
realidad  es  dotar  de  mayor  margen  dinámico  a  los  ND 
más  frecuentes  en  la  imagen,  que  suelen  ser  los  de 
gama  media,  en  detrimento  de  los  muy  bajos  o  muy 
altos,  a  los  que  se  les  roba  parte  de  su  margen  diná¬ 
mico  inicial.  El  resultado  maximiza  el  contraste  de  la 
imagen  sin  perder  información  de  tipo  estructural,  es 
decir,  que  si  en  la  imagen  original  dos  celdas  exhiben 
diferentes  ND,  tras  la  ecualización  seguirán  teniendo 
niveles  de  gris  diferentes,  aunque  con  otros  valores. 

La  función  de  expansión  es  el  propio  histograma  acu¬ 
mulativo  P{z).  Esta  función  se  caracteriza  por  tener 
mayor  pendiente  cuanto  mayor  sea  el  número  de  cel¬ 
das  de  cada  nivel  digita.  Como  se  utiliza  para  reñejar 
en  nivel  digital  original  sobre  el  eje  de  ordenadas  au¬ 
xiliar,  la  separación  entre  los  nuevos  niveles  digitales 
será  mayor  cuanto  los  ND  originales  sean  más  fre¬ 
cuentes. 

z=P{z) 

Una  propiedad  de  esta  modalidad  de  ecualización  es 
que  una  segunda  ecualización  del  histograma  previa¬ 
mente  ecualizado  proporciona  la  misma  imagen,  sin 
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variación  del  contraste,  dado  que  la  abscisa  y  la  orde¬ 
nada  de  cada  punto  de  la  nueva  función  P{z’)  son 
iguales  y  convierten  a  esta  función  en  una  recta  de 
pendiente  unidad. 

Expansión  especial 

Cuando  se  necesite  expandir  con  diferente  amplitud 
dinámica  unas  zonas  y  otras  de  la  imagen  se  puede 
recurrir  a  emplear  una  función  que  mantenga  pen¬ 
dientes  locales  seleccionadas  por  el  analista.  Para  ello 
en  lugar  de  una  recta  se  puede  utilizar  una  polilinea 
de  expansión,  con  una  selección  arbitraria  de  los  pun¬ 
tos  de  cambio  de  pendiente.  De  este  modo  se  consi¬ 
guen  distintos  márgenes  dinámicos  para  las  zonas  de 
interés. 

Binarización 


mica  mediante  una  temperatura  de  corte,  o  en  el  in¬ 
frarrojo  próximo,  mediante  un  nivel  de  brillo  bajo 
que  sólo  aislé  el  agua.  Para  emplear  la  máscara  obte¬ 
nida  se  multiplicará  por  la  imagen  original,  elimi¬ 
nando  de  este  modo  la  zona  no  deseada,  ya  que  las 
celdas  de  esta  zona  en  la  máscara  tendrán  valor  0,  en 
tanto  que  el  resto  valdrá  1. 

Inversión 

La  inversión  del  contraste  se  lleva  a  cabo  mediante 
una  recta  de  expansión  de  pendiente  negativa.  La  in¬ 
versión  puede  realizarse  aprovechando  todo  el  mar¬ 
gen  dinámico  o  se  puede  ajustar  a  cualquiera  de  los 
valores  explicados  anteriormente. 


z 


La  binarización  es  un  caso  extremo  de  la  expansión 
con  saturación  en  la  que  los  valores  m  y  M  se  hacen 
coincidir  en  un  solo  valor  umbral  m  y  se  adopta  una 
función  escalonada.  Desde  este  nivel  hacia  la  iz¬ 
quierda,  todas  las  celdas  adquieren  el  valor  0  y  todas 
las  de  la  derecha,  incluidas  las  que  ostenten  el  valor 
umbral,  adquieren  el  valor  1.  En  la  práctica  la  imagen 
original  queda  convertida  en  una  imagen  binaria  de 
celdas  negras  y  blancas  configurada  a  partir  del  valor 
umbral  elegido. 

,  Jo  V  z  <  M 
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Es  el  procedimiento  para  construir  máscaras  en  una 
imagen.  Si,  por  ejemplo,  se  desea  obtener  una  más¬ 
cara  del  contorno  continental,  eliminando  de  la  ima¬ 
gen  la  zona  marina,  es  posible  construir  una  máscara 
binaria  a  seleccionando  el  umbral  en  la  imagen  tér- 


Expansión  especial  con  saturación  en  los  nive¬ 
les  bajos  y  con  dos  intensidades  diferentes. 
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Binarización  del  histograma. 
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